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ABSTRACT 

In this paper we study correlation functions of circular Wilson loops in higher 
dimensional representations with chiral primary operators of M = 4 super 
Yang-Mills theory. This is done using the recently established relation between 
higher rank Wilson loops in gauge theory and D-branes with electric fluxes in 
supergravity. We verify our results with a matrix model computation, finding 
perfect agreement in both the symmetric and the antisymmetric case. 



August 11, 2006 

*E-mails: sgiombi, rricci, dtrancan@max2.physics.sunysb.edu. 



Contents 



1 Introduction 



2 Kaluza-Klein expansion 



3 Operator product expansion of Wilson loops 



6 
8 

8 

10 

4.1.2 The correlation functionl 11 



4 Brane computation 

4.1 The D3 branel _ 



4.1.1 Coupling to chiral primarie 



4.2 The D5 branel 12 

14 



4.2.1 Coupling to chiral primarie 



4.2.2 The correlation functionl 15 



The correlation functions from the normal matrix model 



16 



5.1 The symmetric easel . 17 



5.2 The antisymmetric easel 18 



6 Conclusion 



22 



A Spherical harmonics and orthogonal polynomials 



23 



1 Introduction 

In the AdS/CFT correspondence local gauge invariant operators are matched with bulk 
supergravity fields evaluated at the boundary of the AdS space pQ[21jS]- An important 
role in the correspondence is also played by non local gauge invariant operators, the most 
notable example being the Wilson loop. 

In M = 4 super Yang-Mills theory Wilson loops are defined (in Euclidean signature) 

as 

W n (C) = ^-^Tr^Pexp jf [iA^r)x» + $/(t)0'|±| ] (1.1) 

where is the gauge field and $/ are the six scalars of the M = 4 multiplet, and 8 1 is a 
constant unit vector in M 6 . The data which characterize the Wilson loop are the shape of 
the integration contour C and the representation 1Z of the gauge group. Supersymmetry 
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restricts C to be a straight line or a circle jlj 1 , while 1Z may be arbitrary. It is a well- 
established part of the AdS/CFT correspondence that Wilson loops in the fundamental 
representation are associated with a classical string surface of minimal area landing on 
the loop jZjJHl- The worldsheet area, being the string infinitely long, is formally infinite 
but the string action can nevertheless be made finite by adding suitable counterterms [H] . 
The expectation value of the Wilson loop is then the partition function of the regularized 
string associated to it. 

Recently circular Wilson loops in representations other than the fundamental have 
been a very active field of investigation. A holographic dictionary has emerged, where 
probe branes in the bulk are related to higher rank Wilson loops in the boundary ^2] [13] 
[Tij |15j 2 . In particular, D3 branes with AdSi x S 2 worldvolume and k units of fundamental 
string charge dissolved on them have been proved to compute expectation values of Wilson 
loops in the rank k symmetric representation 3 . This picture is the natural generalization 
to the AdS$ x 5 5 background of the idea that a fundamental string ending on a D3 brane 
in flat space can be described in terms of a curved D3 brane with a localized spike carrying 
a unit of electric flux, as first proposed by Callan and Maldacena JIB]. On the other hand, 
D5 branes with AdS2 x 5 4 worldvolume and k units of string charge correspond to Wilson 
loops in the rank k antisymmetric representation. Both these branes are half-BPS and 
preserve the same isometries of (jl.lj) . namely 5*0(2, 1) x 50(3) x 50(5). They pinch off 
at the boundary of AdS*, landing on the curve that defines the Wilson loop. 

The intuitive reason for considering these objects is that, to build a Wilson loop in 
the rank k representation, one would start with considering k coincident fundamental 
strings. The D3& and D5& branes can then be thought of as coming from an Emparan- 
Myers polarization effect [Ell [21], which, for k sufficiently large, blows up a 5 2 C AdS$ 
or a 5 4 C 5 5 from the worldsheet of the k coincident strings. This is reminiscent of the 
interpretation of gravitons with large momenta as D3 branes wrapping a 5 3 C 5 5 (giant 
gravitons) or a 5 3 C AdS$ (dual giant gravitons). We can then regard the D3& and D5& 
branes as dual giant and giant Wilson loops, respectively. 

This brane picture has the advantage of automatically encoding the interactions be- 
tween the coincident strings [23 an d yields all non planar contributions to the expectation 
value of the higher rank Wilson loop J2] • 

^dn this paper we consider only half-BPS operators. For loops preserving less supersymmetry see 

2 For 't Hooft loops see and for applications to finite temperature see |17| . 

3 The brane probe approximation breaks down for k much larger than N . In this limit the branes 
backreact deforming the geometry into the supergravity solutions studied in |1L)| and 
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It is well-known that the expectation value of circular Wilson loops in the fundamental 
representation can be computed with a quadratic Hermitian matrix model |22| |23| . It 
has been conjectured that this can be extended to higher rank loops and matrix model 
computations have provided a successful check of the holographic dictionary just discussed 

A small circular Wilson loop, when probed from a distance much larger than its 
characteristic size, can be expanded in a series of local operators of different conformal 
dimension The operators which are allowed to appear in the expansions must pre- 
serve the same symmetries of and therefore must be bosonic, gauge invariant and 
SO(5) invariant. The conformal dimension of some of these operators is not protected by 
the superconformal algebra and therefore they receive large anomalous dimensions and 
decouple in the strong coupling regime. An important class of operators which have pro- 
tected dimensions and appear in the operator product expansion are the chiral primary 
operators. The correlator with a local operator can then be read off from the expansion 
of the Wilson loop j2E! 4 - 

In this paper, we use the D3& and D5& branes to compute the correlation function 
between a circular Wilson loop in a higher representation and a chiral primary operator 
in the fundamental representation. We do this by studying the coupling to the brane 
worldvolume of the supergravity modes dual to the chiral primaries. These modes prop- 
agate from the insertion of the local operator on the boundary to the brane worldvolume 
in the bulk. 

The paper is organized as follows. In sections 2 and 3 we review the operator product 
expansion of the Wilson loop and how to compute correlation functions when the Wilson 
loop is described in terms of a fundamental string worldsheet. To evaluate this one needs 
to study the harmonic expansion on S 5 of the bulk fields which couple to the worldsheet. 

Following the philosophy outlined before, we then replace the fundamental string with 
the D3fc and D5fc branes. We start by investigating the symmetric case in section 4. We 
expand the brane action to linear order in the fluctuations of the bulk fields and find how 
it couples to the relevant supergravity modes. Using the procedure reviewed in section 
3 we compute the correlation function. In the limit of small k we recover the previously 
known result derived using the fundamental string. 

We then move on to the analysis of the antisymmetric case. The D5fc brane now 
extends also in the S 5 directions. Also in this case we compute the correlator between 
the Wilson loop and a chiral primary operator and check that it yields the correct string 
4 For a nice review see |27) . 
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limit. 

As a further check, in section 5, we compare our results against the expressions coming 
from the normal matrix model introduced in this context in [21] and find perfect agreement 
both in the symmetric and antisymmetric case. 

In the appendix we collect some facts about spherical harmonics and orthogonal poly- 
nomials that we have used in the paper. 



2 Kaluza-Klein expansion 

In this section we review the expansion in spherical harmonics for type IIB supergravity on 
AdSz, x S 5 [2E] > and identify the bulk excitations associated to turning on a chiral primary 
operator in the dual Af = 4 gauge theory [29]. These will be later used to construct the 
coupling of the various supergravity modes to the D3& and D5fc branes. 
The Einstein equations read 5 

R = — F - u,F ijkl (2 I) 

yo 

where the 5-form field strength F/g\ is self-dual. In the Poincare patch, the AdS$ x S 5 
solution reads 

ds 2 = — (dz 2 + dx 2 ) + dVL 2 b (2.2) 
-^a i i^2a i 3^4M5 — — 4e^ 1/J2 ^ 3 ^ 4 ^ 5 , F aia2a3a40l5 — — 4e Q , lQ , 2Q , 3Q , 4Q , 5 . (2-3) 

The fluctuations around the background geometry can be parametrized as follows 

Gr m n 9mn h mn (2.4) 

h a /3 = h( a 0) + "g-#a/3, = (2.5) 

V = Kv~y 9 ^ ^U = ° ( 2 - 6 ) 

F = F + 5F, SF ijMm = 5V[ia jk i m] (2.7) 

where h 2 is the trace of the metric on the five-sphere, h 2 = h a pg a ^. Note that the fields 
h uv and h' are related by a d = 5 Weyl shift. To identify the bulk excitation in AdS$ we 

5 In our conventions Latin indices run over the whole 10 dimensional manifold while Greek indices 
(j,, v, . , . and a,(3, . . . run over AdS^ and S 5 respectively. We also choose units in which RAds s = -Rs 5 = 1- 
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expand the fluctuations as follows 6 

h',u = J2 h U x ) Yl (y) (2-8) 

h 2 = YtrtWiv) (2-9) 

-4^e aaia2Q3Q4 fo / (x)V Q y 7 ( 2 /) (2.11) 

where x and y refer to the AdSs and S 15 coordinates respectively, and Y 1 are scalar 
spherical harmonics on the 5-sphere which satisfy 7 

V a V Q y 7 = -A(A + A)Y I . (2.12) 

Spherical harmonics on S 5 can be classified in terms of the SO (6) ~ S , f/(4) R-symmetry 
group. In particular scalar harmonics belong to the [0, A, 0] representation. The fields 
/12 and b appear coupled in the linearized equation of motions. Their equations can be 
diagonalized introducing the linear combinations 

*' = 20(^72) |ft * + 10A,/| (2 - 14) 

which obey the equations of motion 

V M V M s 7 = A(A-4)s 7 (2.15) 
V M V M t 7 = (A + 4)(A + 8)t 7 . (2.16) 

A scalar field in AdS with m 2 = A(A — 4) (with A > 2) transforming in the [0, A, 0] 
representation corresponds to a chiral primary operator Oa of conformal dimension A. 
Therefore, to linear order, the scalar field s 1 corresponds to chiral primaries in the dual 
gauge theory. On the other hand, the scalars t 1 are associated to their descendants, which 
we do not consider in the paper. 

6 We do not consider the harmonic expansion of /i( Q/ 3) as this fluctuation is related to Q 2 Q 2 descendants 
of chiral primaries in the dual super Yang-Mills theory. 

7 We include a brief review of spherical harmonics in the appendix. 
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The linear solutions to the equations of motion turn out to be j2H] 



h 



--Asg^ + 
2 A s g a p 
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(2.17) 
(2.18) 
(2.19) 



A + l 



a 



4f 



V 5 6 



a, 



4 E e < 




(x)v Q y 7 (y) 



(2.20) 



where s = s 7 !^ 7 and 6 = 6 7 y 7 . Using (|2.5|) and the solution (|2.18|) one can identify 
/i 2 = 10 A s. Setting t 1 = in ()2.14jl . one can then deduce s = —b. 

3 Operator product expansion of Wilson loops 

The Wilson loop operator can be expanded in terms of local operators when probed from 
distances much larger than its characteristic size a. For the circular Wilson loop with 
radius a we can write 



In this expression Oi n ) is a local gauge invariant operator with conformal dimension Ar n ), 
and the sum over n runs over both the primary operators and their descendants. This 
operator product expansion must be invariant under the symmetries preserved by the 
Wilson loop. The half-BPS circular loop has ^(r) = 6 1 = const, and therefore preserves 
a 50(5) subgroup of the original 50(6) R-symmetry group. The operators appearing 
in the OPE expansion must therefore contain 50(5) singlets in the 50(6) — > 50(5) 
decomposition. For example, at level A = 2 we can consider the chiral primary operator 
O2 = O^Tr$ 7 $ J , where Cfj is a 50(6) symmetric traceless tensor. Under 50(6) — > 
50(5) it decomposes as 20 — > 1 + 5 + 14 and therefore, containing a singlet, it will 
appear in the OPE of the Wilson operator. A similar analysis can be performed for 
higher dimension operators, which in general will contain covariant derivatives, gauge 
field-strenghts and the fermions of the M = 4 multiplet. Some of them will get large 
anomalous dimension in the strong coupling limit and therefore will decouple. The generic 
expansion looks as follows 




(3.1) 



W(C) 

WW) 



1 + c (2 ) a 2 YX 2) (9)M 2 CfjTr ($ / $ J ) + 
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+c (3 ) a 3 Yf (6)Af 3 C f JK Tr ($'$^*) + C(4) a 3 Tr (6 1 X 1 F + ) +... 

(3.2) 

where (6) are spherical harmonics and J\f n are normalization constants. 

The coefficients appearing in the OPE expansion can be read off from the large distance 
behavior of the two point correlator of the Wilson loop and the local operators 

(W(C)O^(x)) _ a A <") 

(W(C)) ~ C{n) L 2A (") + ' ' ' {6 - 6) 

where it is assumed that the loop radius a is much smaller than the distance L from the 
point of insertion of the local operator. In this paper we will focus only on chiral primary 
operators = C^... /A Tr($ /l . . . $ /a ) 8 . These belong to short representations of the 
superconformal algebra, have protected conformal dimensions, and appear at all orders 
in the expansion ()3.2I) . 

In the AdS/CFT correspondence the chiral primary operators are dual to supergravity 
modes: (9 a corresponds to a scalar of mass m 2 = A (A — 4), which is a combination of the 
trace of the metric and the RR 4- form over S 5 , as we reviewed in the previous section. 

We now briefly discuss the procedure for computing the correlation function of these 
operators with a Wilson loop in the strong coupling regime. The coupling to the string 
worldsheet of the supergravity mode dual to Oa is given by a vertex operator Va, which 
can be determined by expanding the string action to linear order in the fluctuation h^ u 

S = J {fa^detiG^d^dpx") 

= ^7 J d 2 a^J det {g^d a x^d p x v ) + ^{g^d^dpx^h^d^dpx" + 

(3.4) 

The fluctuation of the metric h^ v on AdS^ is given in eq. (j2.17j) . We write the scalar s 1 
in terms of a source Sq located at the boundary 

7 (x, z) = I d 4 x'GA(x'; x, z)sq(x') (3.5) 



s 

where Ga{x';x,z) is the bulk-to-boundary propagator which describes the propagation 
of the supergravity mode from the insertion point x' of the chiral primary operator to 
the point (x, z) on the string worldsheet 

G A (x';x,z)=c ( — ^— — - ) . (3.6) 



3 We take the traces of the chiral primaries in the fundamental representation. 
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The constant c = 2 '2-^n^/ A * s fi xe d by requiring the unit normalization of the 2-point 
function [26J. Since we are probing the Wilson loop from a distance L much larger than 
its radius a we can approximate 

G A (x';x,z)^c^ K , d^^-s 1 , dy ~ A(A ~ 1} /. (3.7) 

The relevant Christoffel symbols are readily computed to be 

1%, = zg^ - -j;St (3.8) 

so that one finally has 

4A 

~ -2A ^ + -^-^s 7 (3.9) 
Inserting this result into (|3.4|) . the coupling to the worldsheet is found to be 

' f dA(-2As)^- = — *— [ dAV A s. (3.10) 



27ra' J 7 a 2 27ra' _ 

In this expression dA is the area element of the classical string. The correlation function is 
obtained from functionally differentiating the previous formula with respect to the source 
s 

{W( -1° I ^ )) = -r' W -4--i 7 /^ <i V^ Ga( f;f,,)*-') 

\yV{C)) 0s {x ) 2na' J 

= -y'(0)_L f dA V A G A (x ;x,z). (3.11) 
Ziia J 

One obtains in the approximations of eq. ()3.7|) 

{W(C)O A (x )) _ nA/2 _ 1 VM a A 

(W(C)) N L2A- [6 - LZ) 

We now move on to studying the operator product expansion of Wilson loops in higher 
dimensional representations. We analyze the rank k symmetric representation first. In 
the bulk this is described by a D3& brane. 



4 Brane computation 
4.1 The D3 brane 

We consider a small circular Wilson loop of radius a placed on the boundary of AdS§. 
The metric on AdS$ can be written in polar coordinates as 

ds 2 AdS = — (dz 2 + dr\ + r 2 # 2 + dr\ + r 2 d0 2 ) . (4.1) 
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The position of the loop is defined by r\ = a and z = r 2 = 0. We take a D3 brane which 
pinches off on this circle as z — > and preserves a SO(2, 1) x 50(3) x 5*0(5) isometry 



The bulk action includes a DBI part and a Wess-Zumino term, which captures the 
coupling of the background Ramond-Ramond field to the brane 

S D 3 = T D3 J v/det(7 + 27ra'F) - T m J P[C W ] (4.2) 

where Td3 = ^ is the tension of the brane, 7 is the induced metric, F the electromagnetic 
field strenght, and P[0(4)] is the pull-back of the 4-form 

C(4) = — j-dri A dip A dr 2 A d</> (4.3) 

to the brane worldvolume. 

We review the brane solution found in [12] . It turns out to be more convenient to use 
a new set of coordinates obtained by transforming {z, r 1; r 2 } into 

a sin rj a cos rj a sinh p sin 9 

cosh p — sinh p cos 9 ' cosh p — sinh p cos # ' cosh p — sinh p cos 9 

In this coordinate system the metric on AdS$ reads 

ds 2 AdS = (drj 2 + cos 2 r] dtp 2 + dp 2 + sinh 2 p (d9 2 + sin 2 9 dcj) 2 )) (4.5) 

sin 77 

where p € [0, 00), 9 G [0,7r], and r\ G [0,7r/2]. The Wilson loop is located at 77 = p = 
0. One can pick a static gauge in which the worldvolume coordinates of the brane are 
identified with {ip , p , 9 , 0} and the brane sits at a fixed point of the S 5 determined by 
the constant unit vector 9 1 G M 6 . The remaining coordinate can be seen as a scalar field, 
rj = rj(p). Because of the symmetries of the problem the electromagnetic field has only 
one component, F^ p (p). In this coordinates the DBI action in ()4.2|) reads 

S DBI = 2N [ dpd9 Sm9s ™ h2 P J cos 2 rj(l+rj'i) + (2m') 2 sin 4 ryF* (4.6) 
J sin rj v ^ 

while the Wess-Zumino term is 

.2 



, cosrysinflsinh p / sinh p - cosh p cos g . 

5vkz = -2iV / dprfg cos 7/ + rj sin 77 — . (4.7) 

sin rj \ cosh p — smh p cos g 

The solution to the equations of motion reads 



1 ikX k\f\ 

sin rj = — smh p , F^ D — ^— , k = —— . (4.t 

' k P ' ^ P 8vriVsinh 2 p' 4A^ V 
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The bulk action has to be complemented with boundary terms for the worldvolume 
scalar rj and for the electric field F^ p [12. . These terms do not change the solution but 
alter the final value of the on-shell action which reads 

S D3 = S DBI + S wz + S hdy = -2N(ny/l + k 2 + sinh" 1 k). (4.9) 

The expectation value of a Wilson loop in the rank k symmetric 9 representation is then 

(W Sk ) = exp (2N(kVI + k 2 + sinh' 1 «)) . (4.10) 

For small k this expression reproduces the result of k fundamental strings 

(W Sk ) ~ e k ^. (4.11) 

4.1.1 Coupling to chiral primaries 

The linearized coupling of the scalar s 1 to the brane worldvolume can be found by expand- 
ing the induced metric on the brane around the AdS$ x S 5 background g mn and keeping 
the first order term in the fluctuation h mn . Since the brane lies completely in AdS§ we 
can write 



Sdbi = T D3 J dV^/det (G^daXW + 2na>F ab ) 



= T D3 J d 4 a^/det(g pi/ d a X»d b X» + 2na'F ab ) ■ 

+ \ {g l u,d a X"d b X v + 2WF afe )~ 1 h pc ,d a X»d b X° + ...) (4.12) 

Here a, b are the brane worldvolume indices. 

The coupling of s 1 to the 4-form in the Wess-Zumino term is obtained by replacing 
C(4) — > C( 4 ) + a ( 4 ) where, using eq. (|2.19jl and the approximation (|H.7)I . the fluctuation 

a (4) is 



so that 



<... M4 - -te^zd's 1 ~ -4A z e^.^s 1 (4.13) 



S$z = ~Td3 f P[a W ) = 4T D3 A J P[C (4) ]s (4.14) 



9 In the limit of N — > oo and A — > oo the symmetric representation coincides with the multiply wound 
Wilson loop |23 ESI 
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where s = ^2 s^V 7 . 

We use now the explicit solution to the equations of motion (|4.8jl to evaluate the 
on-shell value of the fluctuations (j4.12j) and ([4. 14)1 . The first order in the fluctuation in 
(J4~T2|) turns out to be 

9 f , ,/, sin^ / , „ , 1 — sinh 2 p{n~ 2 — sin 2 6) \ 

SZL = AN An 2 / dpdO =- -1 - Ik 2 + ^— — / s. 4.15 

DB/ 7 F sinh 2 pV (coshp-sinhpcosfl) 2 / V 7 

Similarly, the Wess-Zumino term reads 

„»ta a f , sin# / 1 sinh 3 p — sinh p cosh 2 p A „„. 

= 8NAk 4 / dpd# =- 1 + — — " cos# s. 4.16 

7 sinh p \ k cosh p — smh p cos 6 1 / 



The final result for the action is then 
_ c(i) , c(i) 



sinh K sin Q 



S% = S$ BI + S# z = -ANA / dp d9 j— — — a. (4.17) 



o Jo 



(cosh p — sinh p cos 6) 2 



4.1.2 The correlation function 



The prescription for computing the correlation function between the Wilson loop and the 
chiral primary operator is to functionally differentiate the action 1)4.17)1 with respect to 
the source Sq (see eq. (13.5)1 ) 

(W(C)O a (L)) _ 6S& 



(W(C)) 5s 



(4.18) 

s =0 



We approximate the bulk-to-boundary propagator with and use for z the expres- 
sion ()4.4jl . This yields 

(W(C)) " L 2A k a y o P P J Q (coshp -sinh p cos #) 2 + A ' 

(4.19) 

We are neglecting terms of higher order in -j^. 

After performing the two integrals, the final result for the coefficients of the operator 
product expansion turns out to be remarkably simple 

2A/2+1 

Cs feiA = — sinh (A sinh" 1 k) . (4.20) 
v A 



Interestingly enough, this can be expressed in terms of Chebyshev polynomials with ima- 
ginary argument 

(_1)A/2 2 a/2+i J _,y A (m) for A even 
CSkA ~ \ T a (zk) for A odd ( "' 21) 
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where we have used the identities T n (x) = cos(ncos" 1 x) and V n (x) = sin(racos~ 1 x). 

The string limit is recovered when k — > 0. In this regime the S 2 in the brane world- 
volume shrinks to zero size and the D3 reduces effectively to a fundamental string with 
AdS 2 worldsheet. The coefficients ()4.2U|) become 

c SkA z> 2 A / 2+1 VA K = 2 A / 2_1 k (4.22) 
in perfect agreement with the result ()3.12j) found originally in [26J. 



4.2 The D5 brane 

Circular Wilson loops in the rank k antisymmetric representation of the gauge group 
have a bulk description in terms of D5 branes with AdS 2 x S 4 worldvolume and k units of 
fundamental string charge dissolved in them [13] • The D5 description of these Wilson 
loops is valid in the large N, large A limit with k/N fixed. Before moving on to compute 
the coupling of these branes to the scalars s 1 dual to chiral primaries, we briefly review 
the D5 solution [T3| to set up the notation and our conventions. It is convenient to take 
the AdS 5 x S 5 metric as 

ds 2 = cosh 2 u(d( 2 + sinh 2 C# 2 ) + du 2 + sinh 2 w(ch? 2 + sm 2 $d<j) 2 ) + 

+ d6 2 + sin 2 0<ffi| , (4.23) 

where we have written the AdS§ factor as an AdS 2 x 5* 2 fibration. These coordinates are 
related to the usual Poincare patch by 

a cosh u sinh £ a sinh u sin d 

r i = — ; TT^—i — > r 2 



cosh u cosh ( — cos sinh u cosh u cosh ( — cos $ sinh u 

a 



(4.24) 



cosh u cosh £ — cos i? sinh u 

where, as before, a denotes the radius of the Wilson loop. In these coordinates, the Wilson 
loop is at ( — > oo, u = and it is parametrized by ip. The self dual 4-form potential can 
be taken to be 

C (4 ) = 4 ( - - — sinh 4m ] dH 2 A dVt 2 - ( - 9 - sin 29 + - sin 4^ j dtt 4 , (4.25) 

where dH 2 denotes the volume element of the AdS 2 part of the metric. 

Since we want to construct a D5 brane with AdS 2 x 5 4 worldvolume, it is natural to 
take a static gauge in which (, if) and the coordinates of the S A C S 5 are the worldvolume 
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coordinates. Furthermore we can take the following ansatz which preserves the 5*0(2, 1) x 
50(3) x 50(5) symmetry of the Wilson loop 

u = , 6 = const. (4.26) 

and only the component of the worldvolume gauge field is turned on. With this ansatz 
the DBI and Wess-Zumino parts of the D5 action reduce to 



2N 
3^ 



Sdbi = T m J d 6 cr^/det(-f ab + 2ira'F ab ) 

VA J rfCsinhCsin 4 ^l + ^-^^, (4.27) 

S wz = -2ira'iT m J F A P[C {i) ) 

= ^ J dC F^ {^6 - sin 29 + i sin 40^ (4.28) 

where we have used Td 5 = N^/X/Sn 4 and vol(04) = 87r 2 /3. The equation of motion for 
the electric field states that the conjugate momentum is a constant equal to the number 
of fundamental string charge k dissolved in the D5 brane 

„ -i 5C 2N £sin 4 fl 2N /3 „ n 1 . A \ , 

n = — = + [-9- sin 20 + - sin 45 = k 4.29 

2ti5F^ 3tt y/i -E 2 3vr V2 8 J V 1 

where for convenience we have defined E = This equation allows to determine 

the angle 9 at which the D5 sits function of k 

k 

9 k - sin 6 k cos6 k = tt — (4.30) 

while the electric field is given by E = cos#fc. One can check that with this ansatz the 
equation of motion for u is also satisfied. Adding the appropriate boundary terms for the 
electric field and the worldvolume scalars (see ^H] [HD] for details) the on-shell action for 
the D5 brane becomes 

2N 

Sd5 = S DB i + S wz + S hdy = - — vA sin 3 6 k (4.31) 

so the expectation value of the Wilson loop in the rank k antisymmetric representation is 
given by 

(W Ak ) = exp (^V\sm 3 6 k ^ . (4.32) 
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As previously noted in the literature, this result is consistent with the duality between 
the rank k and rank N — k antisymmetric representations: indeed, it can be seen from eq. 
()4.30|) that under k — >• N — k the angle 9 k goes into tt — 9k- It can also be checked that in 
the limit k/N — > 0, in which the 5 4 factor shrink to zero size, (|4.32j) coincides with the 
action of k fundamental strings, as for small k/N eq. ()4.30j) gives 9\ ~ 3nk/2N, so that 
(W Ak ) expky/\. 

4.2.1 Coupling to chiral primaries 

The coupling of the KK scalars s 1 to the D5 worldvolume can be obtained along the same 
lines of the D3 calculation of the previous section. However, besides the fluctuation of the 
AdS§ part of the metric h^, we also need the fluctuation of the metric in the 5 5 direction 
h a p as well as the fluctuation of the 4-form along the 5 4 . The explicit expressions 
can be found in section 2. In particular, in this coordinates the 4-form over the 5 5 is 

a CT1CT2CT3ff4 = 4 sin 4 9 /i(0 4 ) ^ s 1 dgY 1 (4.33) 

where ai,...,a 4 are the coordinates on the 5 4 and /if^) = sin 3 a\ sin 2 <r 2 sin <r 3 is the 
corresponding measure. Differently from the D3 compuation, in this case the 5 5 spherical 
harmonics Y 1 play an active role in the computation since the D5 brane extends into the 
5-sphere. The explicit form of the harmonics is given in the appendix. 

The variation of the DBI part of the action to first order in the fluctuations h^ v and 
h a p reads 

Sdbi = ^f- J Vdet( lab + 2na>F ab ) ( lab + 2ita'F ab y l (h^d a X»d b X» + h a pd a X a d b X?) . 

(4.34) 

Using the explicit solution reviewed in the previous section, it is easy to compute the 
matrix 7 afe + 27ia'F ab . Plugging in the explicit expressions for the fluctuations and using 
the fact that on the D5 solution we have z = a/ cosh£ (this follows from the change of 
coordinate ()4.24|) after setting u = 0), we get after some computations 

S^ BI = nT m [ d(da 1 . . . da 4 fi(n 4 ) sinhCsin 5 ^ ( „ 4A + 8A ) s'Y 1 . (4.35) 

" " J \ cosh ( sin 9 k J 

Performing the integration over the S* 4 , only the 50(5) invariant spherical harmonics are 
selected, namely the harmonics which depends on 9k only, and we get 

S { dbi = ^ [dC smh C sin 5 9 k ( 4A + 8 a) s A Y A >°(9 k ) (4.36) 

on J \ cosh Q sin 9 k J 
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where the suffix on the harmonic indicates that all the quantum numbers except one were 
set to zero by the integration over the 4-sphere. As reviewed in the appendix, these S 4 
invariant harmonics can be explicitely written as 



Y A '°(0 k ) = Af A C%> (cosfl. 



(4.37) 



where (cos9 k ) are Gegenbauer polynomials, and A/a is a normalization constant nec- 
essary to have orthonormality. 

The linear coupling coming from the Wess-Zumino part of the action (|4.28jl can be 
obtained using the expression for the 4-form fluctuation in eq. 1)4.33)1 . and after integrating 
over the S 4 as above, we get 



(i) 
wz 



8NV\ 
3tt 



j rfCsinhCsin 4 9 k cos 9 k s A dg k Y A '°( 



(4.38) 



4.2.2 The correlation function 

The correlator between the rank k antisymmetric Wilson loop and chiral primary opera- 
tor Oa(L) can now be computed plugging ()3.5j) into 1)4.35)1 and ()4.38j) and differentiating 
with respect to the source s A . As before, the bulk-to-boundary propagator can be ap- 
proximated by cz A /L 2A . Recalling that on the D5 solution z = a/ cosh£, the (^-integrals 
can be readily computed and we get 



(W Ak Q A (L)) 
(W Ak ) 



L 2A 



2 A/2 

~3tT 



r AXsm 3 9 k Y Afi (9 k )- 



2 ^VA ( A + l )s , n5 / A cosg 



37rVA(A - 1) 



°' a k) + ^9 ek Y A ' ( 



sin di 



(4.39) 



Using the formula for the derivatives of Gegenbauer polynomials eq. (jA.16)l . we obtain 



AY 



A,0i 



COS0J 



sin 9 1 



de h Y A '°(6 k ) 



Aa 
sin 2 9i 



(AC%\cos9 k ) - (A + 3)cos^C A 2 l 1 (cos^)). 



(2) 



(4.40) 



The correlation function (J4.39)) can then be written as 



(W Ak Q A (L)) 
(W Ak ) 



6(A-2)! 
' (A + 2)! 



2 A/2 

~3tT 



^sin 3 9 k - 



2(A + 1) cosfl fc Cfl^cos^) - AC^'(cos9 k ) 



(2), 



(4.41) 
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where we have factorized out the spherical harmonic evaluated at 6 = 0, Y > (0) = 
A/a^Jt 11 10 - The OPE coefficient c A k ,A we aim to compute is the expression in square 
brackets. Using the recurrence relation eq. (jA.17|) . we find that this expression can be 
written in the compact form 

9A/2 fif A — 2V 

c AkA = — VAAsin 3 fl fc ^ + ^ Ci 2 l 2 (cos 0,)- (4.42) 

This is our final result for the correlation function of rank k antisymmetric Wilson loops 
and chiral primaries. In the next section, we will see that this result exactly matches 
the one obtained from the normal matrix model. As a check, one can verify that this 
expression reduces to the string result of i n the limit k/N — * 0, by using 9\ ~ 37rk/2N 
and eq. (|A.18|) from the appendix. 



5 The correlation functions from the normal matrix 
model 

It is well-known that the expectation value of a circular Wilson loop in the fundamental 
representation of SU(N) can be computed from a quadratic Hermitian matrix model 

(W a ) = ^~J [dM]exp (-^TrM 2 ) ^Tr D e M . (5.1) 

This matrix model is conjectured to capture the physics of the Wilson loop exactly, up to 
instanton corrections [33], to all orders of 1/JV and A. The conjecture extends to higher 
rank Wilson loops as well. The result for the multiply wound Wilson loop has been 
obtained in [12J, whereas [2S] and [H][22] contain the computations for, respectively, the 
symmetric and antisymmetric representations. 

When computing the correlation function between a Wilson loop and a chiral primary 
operator one can substitute the Hermitian model (|5.1jl with a complex one by introducing 
a second matrix M Im and defining z = M + zM Im . In [21] it was shown that, for cer- 
tain representations of the Wilson loop (the multi- winding and the antisymmetric), the 
complex matrix model is equivalent to a normal matrix model, which is a complex model 

10 The OPE coefficient does not include a factor coming from the spherical harmonic evaluated at the 
unit 6- vector 9 1 appearing in eq. Hl.lfl . After a rotation, this vector can always be set to 9 1 = (1, 0, . . . , 0) 
which corresponds to the north pole of S 5 , i.e. 8 = 0. 
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where the matrix is constrained to commute with its conjugate. In the normal matrix 
model the expression for the Wilson loop reads 

For large N, the eigenvalues of this model are distributed in incompressible droplets 
in the complex plane. This leads to interpreting the complex plane as the phase space of 
free fermions, in analogy with the matrix quantum mechanics describing chiral primary 
operators [22] [US]- For example, the Wilson loop in the fundamental representation has 
an eigenvalue distribution given by a circular droplet with constant density 11 





\z\ 




1° 


\z\ 



p{z) = {- ; ' (5.3) 



In it was also shown that the correlation function between a Wilson loop in the 
fundamental representation and a chiral primary operator is given by 12 



(W n O A ) = ^[ [d 2 z] exp (-Tr zz) -j-Trp m Tr z A . (5.4) 

Z N J\z,z\=0 N VAN A 



oA/2 p i 

-/ r 

J N J[z,z]=0 

We now use the normal matrix model to check our results for the coefficients of the 
operator product expansion of higher rank Wilson loops. 



5.1 The symmetric case 

We start by reproducing the result ()4.20|) for c,s k) A using the normal matrix model. Accord- 
ing to the holographic dictionary put forward in [T4"] . we are interested in the correlator 
between a Wilson loop in the rank k symmetric representation and the chiral primary 
operator Oa = ^=^ TrZ . In the limit of large N and large A the symmetric represen- 
tation Wilson loop Ws k effectively coincides with the multiply wound fundamental loop 
, as was shown in [21] and [23]. Therefore we limit ourselves to the simpler case of 
computing (W^Oa), where k is the winding number and corresponds in the brane probe 
picture to the number of fundamental strings dissolved in the brane. 

11 Projecting (|5.3|) into the real axis one recovers the Wigner semi-circle distribution. 

12 The factor 2 A / 2 instead of the 2 _A / 2 of [23] is set in order to have normalizations consistent with 

ESI 
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We start from eq. (4.7) of [21], where we replace everywhere A — > k 2 X 

N 



= r^^ mP i WAV 

X D ' kVAX J 2m y 2Nw J 



(5.5) 



The large winding limit consists in taking N — ► oo while keeping k = ^5 fixed. In this 
limit the integral can be evaluated around the saddle point of the terms proportional to 
N and k 

a, [^w + N\og(l + ^\ i (O.G) 



2iVu; 



which yields 



io* = Vl + k 2 — k. (5.7) 

Inserting to* in (|5.5|) and using 



Vl + k 2 + k = exp(sinh 1 k), Vl + — K = exp(— sinh 1 n) (5-8) 

it is easy to see that 

(W^O A ) = 2N ~ 2 sinh (A sinh" 1 K ) e 2Ar (-^ T ^+ sinh ~ 1 «) . (5.9) 



To get a properly normalized expression one still needs to divide (J5.9J) by 

(W^ ] ) = — *— gaJVCM/T+^+sinh- 1 !^ ( 510 ) 

The final result coincides with eq. (|4.20|) . which we obtained from the brane picture. 
5.2 The antisymmetric case 

To compute the OPE coefficients of Wilson loops in the rank k antisymmetric represen- 
tation, we have to evaluate the following correlator in the normal matrix model 



oA/2 k\/8N f 1 __ 

{w - °*> = i^vs L^«;^ w TrA (5 ' li; 



Z N N A / 2 VA J [z ,- z]=0 

This matrix integral can actually be solved exactly, as was shown in [21] , and similarly to 
the case of the fundamental representation, it can be written as a fc-dimensional contour 
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integral. However, it does not seem to be easy to take the large N and large k limit with 
k/N fixed from such an expression. Here we follow a different approach to get the above 
correlator in this limit. First, as in we find it convenient to rewrite the trace in the 
antisymmetric representation using the corresponding generating function 



Tr 



±(z+z) 



dt 
2m 



t^exp 



Tr log 1 



1 i 
- e2 
t 



(5.12) 



Since we expect the correlator to be real, it is also convenient to replace Trz A — > ^(Trz A + 
Trz A ). The idea is then to view the insertion of the chiral primary Trz A in (J5.11|) as a 
small perturbation of the gaussian potential, by writing 



[d 2 z]e-^ z ~ z) Trz A e TAo ^ e ^ W 



(z+z) 



z,z]=0 



= Z N 

= z N 







1 



da \Z N (a) ./ M= o 
d 



^2 z j e -Tt(«g)+f ( Tr2 A +T r z A ) e T r log(l+i V 7T^+ Z >) 



da 



exp 



Trlog ( 1 + ] - e ^^ {z+z) 



a=0 

(5.13) 



a=0 



where we have introduced an a-dependent partition function 



Z N (a) 



2 ] -Tr(zz)+| (Tr2 A +Tr2 A ) 



[z,z]=0 



[d 2 z}e 



(5.14) 



and we have used that Z^(a) = Zn + 0(a 2 ) 13 . The problem is now to evaluate the 
correlation function (J5.13)) in the normal matrix model with the deformed potential 



V(z, z) = -Trzz + |Tr(^ A + z A ) 



(5.15) 



Normal models with potentials of this kind were previously studied in the literature, for 
a recent account see for example |S3]jSS]. To solve the model at large N, one can as usual 
go to the eigenvalue basis at the expenses of introducing a Vandermonde factor, and 
determine the eigenvalue density p a (z, z) in the continuum limit. The density is found by 
solving the saddle point equation 14 



z z 

2 



(5.16) 



13 



This follows from the fact that in the matrix model with gaussian potential (Trz z 



(Trz 1 



0. 



14 As in |2S|, the term exp Trlog ^1 + \ e = 
leading order at large N. 



does not modify the saddle point equation at 
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where the term in the right hand side comes from the Vandermonde factor. Once the 
density is known, the correlation function in ()5.13j) becomes 



cxp 



Trlog ( 1 + - ea 



exp 



1 i 

1 + -e% 



{z+z) 



(5.17) 



It is known that for potentials of the kind V(z, z) = —zz + f(z) + f(z), the density is a 
constant (equal to in the normalizations we are using here) inside a certain droplet in 
the complex plane and zero outside. For the gaussian potential, as reviewed previously, 
the droplet is just a circle of radius y/~N (to compare with eq. ([5.3)1 . one has to rescale 
:z), while the term proportional to a induces a deformation of the circle which 



preserves its total area (since we do not change the number of eigenvalues). It is not 
difficult to find the shape of the droplet which solves the saddle point equation (j5.16|) at 
leading order in a. It is convenient to work in polar coordinates z = re 1 ^ . The curve 
which bounds the droplet can then be written at first order as 



r(0) = VN(l + a /(<!>)). 
Clearly /(0) has to be periodic, and may be written as 

oo 

f{4>) = 5Z a ™ cosn< l ) 



(5.18) 



(5.19) 



n=l 



where only cosines appear because of the symmetry of the potential (|5.15|) under z <-» z, 
and the mode with n = is excluded by requiring the area to be preserved. The saddle 
point equation now reads 



re - 



A« A-l -t(A-l)* = 1 

2 7T 



2k 



/]V(l+« /(<//)) 



dr' ■ 



r e 



(5.20) 



Expanding the integral at first order in a and plugging in the Fourier expansion (|5.19J) . 
we see that this equation is satisfied if a& = iV A//2_1 ^ and all other a n vanish, so we find 
that the shape of the deformed droplet is given by the curve 



(5.21) 



(0) = Vn (l + | A N^ 2 - 1 cos A0) . 



Before moving on to compute (|5.13|) . we can check the validity of the method by 
applying it to the computation of the correlator (|5.4j) when the Wilson loop is in the 
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fundamental representation. In this case, following the same steps as above, in the large 
N limit we arrive at 

2 A/2 q 



(W a O A ) 



AiV A / 2 da 
2 a/ 2 Q l 



d 2 zp a {z,z)e^^~ z) 



a=0 



da Nn 
2 A / 2 V / A 1 r2n 



2k 



/]V(l+f A jV A / 2_1 cos A</>) 



drreVx rcos4> 



N 2vr J 



oA/2 fX 

^ cos *cos A<f) = * Ja(VA) 



a=0 

(5.22) 



which is the result first found in |36| and the correct large N limit of the exact formula 
(B (with = 1). 

Going back to the antisymmetric representation, we have to evaluate 



dt . , d 
—t* '^-exp 
27T2 eta 



N d z p a (z, z) log ( 1 H — e 



1 



a=0 



2ni da 



X j tl 1 (.:..:) log ( 1 + - e^V 7 !^) 



X 



Q = 



x exp 



iV / d 2 z p (z,z)\og ( 1 + - e 2V iv 



(5.23) 



where in the last line po is just the circular droplet density. Since the exponent is inde- 
pendent of a, the t integral can be evaluated in the supergravity limit of large A exactly 
as in we first make a change of variables t = e^ w , then the saddle point of the 
exponent is found to be 



cos6i 



(5.24) 



where 9^ is defined as in eq. f!4.30|) . The exponent in ()5.23|) gives a term proportional to 
the expectation value of the Wilson loop, while the prefactor is evaluated at the saddle 
point. After dividing by (WA k ), the OPE coefficient can then be obtained as 

(W Ak A ) 2 a / 2 d 



(W Ak ) VAN " A / 2 da 
2 A / 2 iV v / A d 2 



AiV A/2 da 7T J 



N j d 2 z p a {z, z) lo. 

1+f AiV^-icosA^ 

cos Ofc 



i + e ^m~ 2c ° s9 ^ 



drr(r cosd> — cos 9i 



a=0 



a=0 

(5.25) 



where the lower limit in the r integral comes from the fact that in the large A limit the 

integral has support only in the region rcos0 > cos^ 15 . After doing the derivative, 

15 The upper limit in the <fi integral is rigorously 9k + 0(a), but it is easy to see that the correction 
does not contribute at first order in a. 
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(|5.25jl gives the final result 



(W Ak Q A ) = 2 A / 2 VA\ 
(W Ak ) 




I dcj) cos Acf) (cos0 — cos^). 





(5.26) 



Remarkably, the integral in (|5.2(jj) precisely reproduces the Gegenbauer polynomials ari- 
sing in the bulk computation, and the final result is 



6 Conclusion 

In this paper we computed the correlation function between a higher rank Wilson loop 
and a chiral primary operator in the fundamental representation using branes with electric 
fluxes. Following the proposal of we considered a D3fc brane for the rank k 

symmetric case and a D5& brane for the antisymmetric one. We then checked our results 
with the normal matrix model discussed in j21j finding perfect agreement in both cases. 

We focussed on chiral primary operators but it should not be difficult to extend our 
computation to operators corresponding to other supergravity modes. For example, the 
KK modes of the dilaton are necessary to compute correlation functions of Wilson loops 
and Tr$ A F2. 

It would be worthwhile to study more general representations of both the Wilson loop 
and the chiral primary operator. A particularly interesting issue to address would be 
understanding from our brane picture the selection rule found in [21]: for Wilson loops in 
the rank k antisymmetric representation the only non vanishing correlators involve chiral 
primaries with traces over Young diagrams with at most k hooks. Another direction to 
pursue may be considering the correlation function between higher dimensional Wilson 
loop and a chiral primary operator with A ~ N. In the bulk this would require to 
study the bulk-to-bulk exchange of supergravity degrees of freedom between the electric 
branes describing the Wilson loop and the (dual) giant gravitons associated with the chiral 
primary. 
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A Spherical harmonics and orthogonal polynomials 

In this appendix we collect some facts about spherical harmonics and orthogonal polyno- 
mials we have used in the paper. We follow the treatments of jHZ| and [38J. 

Spherical harmonics in d dimensions are eigenfunctions of the Laplacian on the unit 
(i-sphere 

V^ J (0) = AK'(n) (A.l) 

where the Laplacian is 

V 2 (d) = — L= dt g* dj (A.2) 

with the metric given by gij = diag(l, sin 2 0^(1, sin 2 0d-i(- ■ •)))• The integer multi-index 
I = (Id, ... ,h) satisfies 

Id > ld-i > ■■■ > h > \h\. (A.3) 
The general solution to eq. (|A.1|) is 

Y l ^(e d ,...,e 1 ) = e ^f[nP!:- 1 (e n ) (A.4) 

v n=2 

where we defined 

n P l L {9) = n ci(Em9)-^^P£^% m (coBe). (A.5) 
In this expression P~ m (x) is the Legendre function of the first kind and the constant 

' (2L + n - 1) (L + I + n - 2)P 1/2 



(A.6) 



2{L-l)\ 

is chosen to ensure the ort honor malizat ion condition 

fi(n d ) Y I Y 1 ' = 5 11 ' (A.7) 



where ix(VL d ) is the measure over S d . The integration over S d 1 selects only SO(d) invariant 
harmonics 

" ' " ' (A.s) 



J I Id 
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The eigenvalue A depends only on l d = A because of the 0(d + 1) symmetry of the 
problem and it can be found by studying the action of the Laplacian on SO(d— 1) invariant 
spherical harmonics 



d 



sin^" 1 6 d 86 d 



OVdJ 



(A.9) 



After the change of variable x = cos 9 d: this is recognized to be the Gegenbauer equation 

d 2 



(1-x 2 ) 



dx 2 



dX Tx) Y A '°(-) = ^Y^(x). 



(A.10) 



The solution to this equation is 

A A = -A(A + d-l), 



Y^°(x)=M A C K A 2 >{x) 



(A.ll) 



( ) 

where C ' A 2 are Gegenbauer polynomials and the constant A/a can be obtained from the 
orthonormality of the F A,0 's 



A/a 



A!(2A + d-i) [r(^)] 2 r(f) 



2 4 - d 7r^r(A + d-i) 



1/2 



(A.12) 



The Gegenbauer polynomials (x) are a generalization of the Legendre polynomials 
and can be obtained from the following generating function 



1 y 
(l-2xt + tr=^ C - {x)tA - 



(A.13) 



We list the first few of them 



C& x \x) = 1 

C {X \x) = 2Xx 

C (x \x) = -A + 2A(1 + X)x 2 

C (x \x) = -2A(1 + X)x + ^A(l + A)(2 + A)x 3 . 



They satisfy the normalization condition 



J' dx(l - x 2 f- 1 ' 2 



C 



,(A) 



I 7T 



T(A + 2A) 



(A + A)r 2 (A)r(A + i) 



(A.14) 



(A.15) 



for A > -1/2. 
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In the paper we have used the following formula for the derivative of a Gegenbauer 
polynomial 

(1 - x 2 )d x C [ £\x) = -AxC<£\x) + (A + 2A 
and the following recurrence relation 

AC^\x) = 2(A + A - l)xC^_ x {x) - (A + 2A 
We have also used that 

° A A!(d-2)! • 
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